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Abstract 

We study the dynamical evolution of a system with a phase space consisting 
of configurations with random energies. The dynamics we use is of Glauber 
type. It allows for some dynamical evolution ang aging even at very low 
temperatures, through the search of configurations with lower energies. 
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Lots of efforts are being devoted to the study of the out of equilibrium dynamics of 
spin-glasses, in order to understand several experimental findings, like aging, and the joint 
existence of dynamics and memory in temperature cycling experiments JT|. Various ap- 
proaches have been used, including some phenomenological studies based on droplet or 
domain growth 0, numerical simulations |3|, mean- field models j§-§|, diffusion on tree like 
structures @, and models based on the existence of traps in phase space. 

In this note we shall elaborate along the lines of this last approach. It deals directly with 
the structure of the phase space, made of many metastables states, figuring the configurations 
of a spin system. Such a picture has been put forward in ||, and generalized in |J, presenting 
the phase space as a random energy landscape made of traps, with a broad distribution of 



trapping times: the energies are the low lying ones of a REM ||10|| , so they have an exponential 
probability distribution, and the trapping times, given by an Arrhenius law, have a power- 
law distribution with infinite mean. The dimension of the space is infinite (equivalently, all 
the traps are connected). In this frame, the diffusion is anomalous PH|IT|| , and aging is 
present. One virtue of this approach has been to point out a simple "kinematic" ingredient 
which induces aging: in a system with a broad distribution of trapping times, the probability 
to be at time £ in a trap of lifetime r always depends on t, and this induces aging. Besides 
this kinematic effect, there might well exist another, "dynamical", source of aging, namely 
the explicit evolution with time of the distribution of trapping times. We shall provide 
hereafter such an example of a dynamical aging process. 

In the trap model, the probability of hopping from one configuration % to another con- 
figuration j, Wi^j, depends only on the energy of the i: the energies are in fact seen as 
energy barriers, which are then uncorrelated. Another case was studied in ||12|| , where Wi^j 
depends only of the arrival configuration j, which allows a great simplification of the mas- 
ter equation. Only these two extreme cases have been studied so far. Hereafter we shall 
consider a case where the transition probability, which is of Glauber type, depends on both 
configurations, so the energy barriers are now correlated: this seems a priori a more refined 
way to define a dynamic for the REM. Furthermore this case allows for the existence of a 



dynamical evolution even at zero temperature, through the search of configurations with 
lower energies. During this evolution it becomes more and more difficult for the system to 
find a lower configuration, which results in a slow down of the dynamics. We shall show 
through an explicit solution at zero temperature that this mechanism gives rise to an aging 
effect in which the system never reaches equilibrium. This aging does not take its origin in 
energetic barriers, but rather in some kind of entropic barrier, namely the low probability 
of finding a favourable direction in phase space. Another example of aging due to entropic 
barrier has been proposed recently by Ritort [|l3fl . 

The model we consider is defined as follows: The system can be in any of N configurations 
i = 1,..,N. The configuration % has an energy E^. The energies are independent random 
variables with distribution P(E). The probability of hopping from one configuration to 
another can be defined in several ways, the only a priori constraint being the detailed 
balance: 

Wi^e-P^ = W^e-^ (1) 

For example, for the traps model, where the lifetime of configuration % is Tj = exp(—j3Ei), 
the transition rates are 

w _ = eg^_ Here 

we consider a transition probability depending 
on both Ei and Ej, given by the Glauber dynamics: 

Wir -' = N 1 + expiPiEj - Ei)) (2) 

As mentioned before, this system is quite different from the trap model (see figure |l|). For 
instance, at zero temperature, the jump to a lower state is allowed in our model, while it is 
impossible to jump out of a trap. 

In this study, we will be interested in computing the law of diffusion (the number of 
configurations reached at time t), the evolution of the mean energy with time, the probability, 
given a time t, to be in a configuration of lifetime r (we shall define the lifetime r precisely 
below), which we will note Pt(r), and the two-times correlation function C(t w + t,t w ). This 
last quantity is defined as the mean overlap between the positions of the system at times 



t w and t w + t: the overlap is simply either 1 if the system is in the same configuration, or 
if it has moved. Assuming an exponential decay out of the configurations, the correlation 
function is related to ptij) through: 

C(t w + t,t w ) = / dr Vtw {T)e- ttT (3) 

J 

Before turning to the exact solution of the master equation at zero temperature, it is 
useful to start with a discussion of the trapping time distributions. When the system is in 
configuration i, with energy the probability of going away per unit time is 

Ps (E i )=Y, W i->3 ■ ( 4 ) 
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For large N and using the definition (0) of the transition probability, one gets at zero tem- 
perature: 

p a (Ei) = f El dE' P(E') (5) 



The "trapping time" Tj is defined as } w , . It depends only on the energy of the configuration, 
through the relation: 

We deduce that, regardless of P(E), the a priori distribution of the lifetimes is: 

dr 

P (r)dr = —6(r - 1) (7) 

As in the trap model, this is a broad distribution with a divergent mean lifetime (al- 
though here it is just marginally divergent, for any P(E)). As was shown in |§, this fact in 
itself creates an aging effect. In our case there is an additional effect because the effective 
distribution of lifetimes evolves with time. After k jumps the system will be in a lower en- 
ergy configuration, and the probability Pk{r) of having a lifetime r is different from Pq{t). 
The zero temperature dynamics gives the recursion relation: 

P k+ i(r) = Po{r) J dr' r' P k {r')6{r - r') (8) 



which leads to: 



W = 2^V-U. (») 



The typical lifetime thus increases exponentially with k (it means that the diffusion is 
logarithmic). This will add up to the usual effect of a diverging mean lifetime in order to 
induce aging. 

We now proceed to the solution of the master equation at zero temperature using the 
Laplace transform. Denoting by Pi(t) the probability of being on configuration % at time t, 
we have 

= E T iJ Pi® (io) 

with Tij = Wj_>i for i ^ j, and J2i Ty = 0. The Laplace transform pi(4>) = J dtpi{f)e~^ 1 
satisfies the equation: 

where the lifetime r% is defined as before by ^ = Tji and where we will take Pi(0) = jj. 
To solve this equation we introduce the Laplace transform of the occupation probabilities 
for all configurations of energy E: 

ffrfl^pMWE-E,) . (12) 
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Using ([TT| ) one derives for / the equation: 

f(E, 0) = g(E, 0) (l + jH dE'f(E', 0)) (13) 



where 



_ , N 1 _ S(E — Ej) P(E) , tS 



j r 1 tj V 1 r(S) 

The self consistency equation ([H|) is easily solved and gives: 



/(£, 0) = <?(£, 0) exp (jT g(E', 0)d£') . (15) 



We can now use this solution to compute the physical quantities of interest. We start 
with the probability Pt(r) to be a t time t in a configuration of lifetime r. Its Laplace 
transform with respect to t is given by: 

Mr) = Xrtmr - x.) = j—^ s(E(Tl fl = IrT ^, (16) 

This Laplace transform can be inverted and gives: 

ft (r) = tT ~ r l + T exp (-1) 0(r - 1) (17) 

We see that this expression decreases as for t « r, as for a model of traps for Po( r ) = 
1/t 2 , but the exponential term makes the probability of being at time t in a configuration 
with lifetime smaller than t very small. 

The correlation function can also be computed for large t and t w : 

C(t w + t,t w )~-^- (18) 

We see immediately the j- scaling of the correlation function. The behaviour lim^oo C(t w + 
t,t w ) = is a consequence of the existence of dynamics even at zero temperature, while the 
behaviour in the other limit, lim^^oo C{t w + t,t w ) = 1, reflects the existence of "weak 
ergodicity breaking" (the value 1 for the limit is a consequence of T = 0). 

Let us emphasize that all these results are independent of the distribution P{E). The 
two main hypotheses of the derivation are the fact that the connectivity is infinite, and the 
temperature has been taken equal to zero. These results can be partially extended in the case 
where P(E) is an exponential distribution, P(E) ~ pexp(pE)9(—E). Such a distribution 
has been found in mean field spin glass models |14| , and it is at the heart of the trap model 
description, since it leads to a broad distribution of lifetimes when (3 > p. Specializing to 
this case, we can first explain in more details the zero temperature dynamics studied above. 
Indeed, the relation between energy and trapping time can be explicited as: T{ = e~ pEl and 
the energy distribution evolves then as 

Pk(E) = i^L exp( P E)6(-E) (19) 
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The mean energy decreases as —k/p with the number of visited configurations, or equiva- 
lently as — logt/p (remember that the diffusion is logarithmic). 

For finite temperature, the set of self-consistent equations is more complicated; we write 

dj = exp(j3Ei), so that 

±=«ff'"^ (20) 



T { 1 JO 1 + vP/P 



and we obtain 



PiW = J^T ( 21 ) 

poo 

/(A, 0) = g(X, 0) + / dp g{\ + fi, 0)/(p, 0) (22) 
Jo 

Taking A = r®l p and writing g(r,(f)), /(r, 0) instead of g(r^' p ,(j)), /(r^/ p ,0), we obtain 
for / the following scaling: 

/(r,0) = ^r 1 -f/ i (0r) (24) 



with /i(a;) behaving as -\ for a; >> 1 and J °° dxh(x) finite. After some calculations, it can 
be shown that p t (r)dr behaves like ^ for 1 << t « r, and as ^ ( ~) " for 1 << r << t, 
and for (3 » 1. 

We obtain thus qualitatively the same behaviour for p t {r)dr as before, and also for the 
correlation function: the dynamics is not modified by a small temperature. 

Note that this behaviour holds in the limit of infinite N\ for any finite N the system 
eventually thermalises, after a time proportional to N (for example, the minimal energy of 
iV states with exponentially distributed energies is —l^sK so ^ takes a time N to find it) 

For this model, some of our results are similar to those of ||: the mean energy decreases 
as —log(t), and at time t the most probable configurations are of lifetime r = t. Nevertheless 
we must emphasize that the mecanism is totally different: in a model of traps, the mean 
energy decreases because the system visits more and more traps, and so it has more and 



more chances to find deep ones. At each step Pk{E) remains the same: = Pq. The 
diffusion is in t p/2/3 , and the energy at step k evolves like —log(k), because the minimum 
of k energies distributed according to the exponential distributionis in —log(k). In our 
model, on the opposite, the energy distribution that the particle sees evolves, and so does 
the distribution of trapping times (see ([5])); the energy decreases in fact because it is easier 
to find a configuration with lower energy (it takes a time proportionnal to exp(-pE)) than 
to move by thermal activation (a time exp(—f3E) is needed). It means that the system 
spends less time in a given configuration, but, since the energy at step k decreases like 
—k, the diffusion is much slower (logarithmic instead of a power-law), so we finally get the 
same behaviour for the mean energy as a function of time. However the main feature is 
that there exists a zero-temperature dynamic, which is qualitatively not modified by a small 
temperature. 

It would of course be very interesting to be able to generalize this approach beyond the 
case of an infinite connectivity. For instance a more realistic definition of the REM dynamics 
could be to start from the definition of the REM in terms of spins with p (— > oo) spin 
interaction |TU|,|rSJ, and use the transition matrix resulting from single spin flip dynamics. 
This seems rather complicated at the moment. 

We thank J.P Bouchaud, R. Burioni, L. Cugliandolo, D. Dean C. De Dominicis and J. 
Kurchan for many helpful discussions on related topics. 
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FIGURES 



FIG. 1. (a):traps model; (b): "steps" model; remember that the connectivity is infinite, so that 
such a picture can be misleading! 
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